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Abstract 

Projection Monte Carlo calculations of lattice Chiral Effective Field Theory suffer from sign 
oscillations to a varying degree dependent on the number of protons and neutrons. Hence, such 
studies have hitherto been concentrated on nuclei with equal numbers of protons and neutrons, 
and especially on the alpha nuclei where the sign oscillations are smallest. Here, we introduce the 
“symmetry-sign extrapolation” method, which allows us to use the approximate Wigner SU(4) sym¬ 
metry of the nuclear interaction to systematically extend the Projection Monte Carlo calculations 
to nuclear systems where the sign problem is severe. We benchmark this method by calculating 
the ground-state energies of the ^^C, ®He and ®Be nuclei, and discuss its potential for studies of 
neutron-rich halo nuclei and asymmetric nuclear matter. 

PACS numbers: 21.10.Dr, 21.30.-x, 21.60.De 
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I. INTRODUCTION 


Lattice Chiral Effective Field Theory (EFT) is an ab initio framework [l-5| which has 
recently been applied to studies of the structure of ligM and medium-mass nuclei 6|-l9|, 


as well as to the physics of dilute neutron matter 


Hln|. 


In particular, t 


le structure of 


12 ( 


C and has recently been elucidated using lattice Chiral EFT 


12 


141. The role of 


al ong with its anthropic implications for 


, |l6[. These successes notwithstanding, 


the Hoyle state in ^^C has also been investigated, 
the viability of carbon-based life as we know it ll5 
lattice Chiral EFT has so far mainly been applied to “alpha nuclei”, i.e. to nuclei with A 
a multiple of 4 and with an equal number of protons and neutrons. This limitation is due 
to the appearance of sign oscillations (the so-called sign problem) in the Projection Monte 
Carlo (PMC) calculation at large Euclidean time. The sign problem also necessitates the use 
of relatively large lattice spacings of a ~ 2 fm, in order to moderate the repulsive short-range 
contributions in the leading-order (LO) lattice Chiral EFT Hamiltonian. Even then, the low- 
lying spectra of most alpha nuclei can only be extracted after considerable extrapolation in 
Euclidean time. For instance, this has so far precluded studies of neutron-rich halo nuclei 
where sign oscillations are more severe. 

In spite of the sign problem prevalent in lattice Chiral EFT, useful results have been made 
possible by the observation that nuclei can be approximately described by a Hamiltonian 


which respects t 
interchangeable 


;he Wigner SU(4) symmetry where spin and isospin degrees of freedom are 
Since Euclidean time projection with an SU(4) symmetric Hamiltonian 
is possible for most nuclei without a sign problem, one can obtain a trial wave function which 
is much closer to the ground state of the full Hamiltonian. The shorter Euclidean projection 
time possible with the full Hamiltonian then becomes sufficient, and moreover the coupling 
constant of the SU(4) symmetric Hamiltonian can be varied in order to generate a large 
set of independent trial states. Such self-consistent extrapolation in Euclidean time, also 
referred to as “triangulation”, allows for a much more precise determination of the properties 
of the nuclei under investigation than would otherwise be possible. The usefulness of the 
Wigner SU(4) symmetry was also noted earlier in a different but related context. It was 
shown in Ref. ISj , that the nuclear interactions are SU(4) symmetric in the limit of infinite 
colors, Nf^ —)■ CX3. Further insight was obtained in Ref. 19|, where it was found that in the 
limit of large R-wave scattering lengths, the two-nucleon interactions exhibit Wigner SU(4) 
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symmetry. 

We shall now take the idea of the SU(4) symmetric auxiliary Hamiltonian one step further 
by considering a weighted sum of the full and SU(4) symmetric Hamiltonians. In particular, 
the weight of each component in the Hamiltonian is controlled by a parameter dj^ such 
that for = 1 we have the full Hamiltonian only, and for = 0 we have the SU(4) 

symmetric Hamiltonian only. For values 0 < d^^ < 1, we have a linear combination of the 
two Hamiltonians. The properties of the physical system are then found by extrapolation 
to d^j —)■ 1, combined with an extrapolation in Euclidean time. As will be shown here, this 
allows us to circumvent the sign problem to a large extent, and opens up the possibility 
to study neutron-rich nuclei and systems where N ^ Z. It also allows us to access much 
larger Euclidean projection times with the physical Hamiltonian, thereby greatly increasing 
the level of conhdence in the Euclidean time extrapolation. We shall demonstrate this 
“symmetry-sign extrapolation” (SSE) method by refining our earlier results for the ground 
state of and by computing the ground state energies of ®He and ®Be. 

We begin in Section [TTl by introducing the SSE method and discussing the origin of the 
sign problem in lattice Chiral EFT, and in Section HTTl we present our method of analyzing 
the PMC data, along with our updated results for ^^C. In Section HVl we apply symmetry- 
sign extrapolation to the ®He nucleus and its mirror isobar ®Be. Finally, in Section IVl we 
discuss the prospects for extending our method to more neutron-rich nuclei and to nuclei 
where N ^ Z. 


II. THE SYMMETRY-SIGN EXTRAPOLATION METHOD 


In order to introduce the SSE method, we hrst briefly recall the ingredients of lattice 
Chiral EFT at LO in the EFT expansion in Q/A. The LO contribution is treated non- 
perturbatively, while NLO and higher order terms are included as a perturbative correction. 
We use the same notation as in Ref. j4|, and further details of the lattice action can be found 
there. At LO, the lattice Chiral EFT partition function is given by 


2lo = / exp 1-S„«)1 Tl|A4o«,A',-l)"-A4oW.0)i-. 


( 1 ) 


f-Ei. 


where vrj- is the pion field and is the free pion lattice action. As in Ref. pj , we are using 
an improved LO operator where the contact interactions depend on the momentum transfer 
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through a smooth smearing function /(^. The normal-ordered LO transfer matrix operator 
at Euclidean time step = 0,..., is 

= ■ exp [-H^o{^:'I,n^)a^] (2) 


with 



Hiree + ^ /(^ : Ps®Ps(-^ ^ 

q q,S 

+ Icp /(^ : : + \Cs2j2 Y • Ps,i{^Psj{-^ ■ 

qj q,S,I 


9a 


V A 


s'^i 


t) Ps. 


n 


(3) 


where the spin vector S and the isospin vector / indices range from 1 to 3, the parameter 
Of = a^/a is the ratio of temporal and spatial lattice spacings, is the lattice gradient 
along direction S', and = a^{m^ + Q). The operator p is the total nucleon density pg 

is the spin density AIct^A, pj is the isospin density AVj-A, and p^j is the isospin density 
N^'agTjN. The couplings C, 0^2, Cj 2 and 0^2 j 2 satisfy 


C=-3Cs2,I2 = ~{Cs2+Cj2), 


(4) 


such that the smeared contact interactions only contribute to even-parity channels, where 
we have antisymmetry in spin and symmetry in isospin (or vice versa). 

The PMC calculations are performed with auxiliary fields coupled to each of the densities 
p, pg, pj, and pgj. The LO auxiliary-held transfer matrix at Euclidean time step is 


^Loln^is, Ss, Sj, Ssj, n'j) = : exp <{ Y n^)psj{n) 




n,S,I 


+ E s{n,n^)p{n)+ i^/CY^tY^sin,n^)ps{n) 

n n,S 


+ i^/Cj2atY^i(^^^t)Pii^) + '^\JCs2j2atY ^sA^^^t)PsA^) f O (5) 

nj n.S.I I 


where the physical values of the LO operator coefficients are such that all factors inside the 
square root symbols in Eq. (I5|) are positive. We now also dehne an SU(4) symmetric transfer 
matrix. 


mY = ■ exp [-H^a^] :, 


( 6 ) 
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with 


H 4 = H^ree + ^ /(^ ^ : , (7) 

where 6*4 < 0 is the associated coupling constant. Again, we can rewrite the interaction 
using auxiliary helds, giving 

^irauU-s) = : exp I -iffreetti + \/^C^ ^ s{n, nt)p{n) | : . (8) 

Let us now consider the projection amplitude we obtain for an A-body Slater determinant 
initial state 


Id') = iV'i) A IV' 2 ) A---A IV’a), 


(9) 


with an auxiliary-held transfer matrix M^ux ■ The projection amplitude is given by det(iW), 
where M is the Ax A matrix we obtain from the single nucleon amplitudes 




( 10 ) 


Let us dehne W*[M] as the set of unitary matrices such that U'^MU = M*. It can be shown 
that det(iVT) is positive semi-dehnite if there exists some antisymmetric matrix U G W*[M]. 
The proof follows from the fact that the spectra of M and M* must coincide, and the 


real spectrum of M must be doubly degenerate as a result of the antisymmetry of U 


20 |. 


A straightforward generalization of this result is that the projection amplitude det(AT) is 
positive semi-dehnite if there exists a unitary operator U G W*[Miux ] for all = 0, • • • , 
and if the action of U on the single-particle states • • •, I'^n) can be represented as an 
antisymmetric Ax A matrix. 

We note that contains the spin and isospin matrices CT 2 and T 2 - PMC calcu¬ 

lations with the SU(4)-symmetric theory are then free from sign oscillations whenever the 
initial single nucleon states are paired into spin singlets or isospin singlets Qfl- We may 
then dehne the “interpolating Hamiltonian” H as 


H — dj^H^Q + (1 


( 11 ) 


which depends on as well as the (unphysical) coupling constant C 4 . This can also be 
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viewed as giving the interaction parameters a linear dependence on 


C{dh) = df^C + (1 — di^)C^, 

Cs2i2{dh) = df^Cg2j2., QAidh) = 


( 12 ) 


By taking d^j < 1, we can always decrease the sign problem to a tolerable level, while simul¬ 
taneously tuning 6*4 to a value favorable for an extrapolation d^^ —>■ 1. Most significantly, 
we can make use of the constraint that the physical result at d^^ = 1 should be independent 
of C 4 . The dependence of calculated matrix elements on df^ is smooth in the vicinity of 
dj^ = 1. In what follows, we shall explore the properties of H for various nuclei of physical 
interest, and determine to what extent it can be used to circumvent the sign problem, which 
at d^ = 1 becomes exponentially severe in the limit of large Euclidean time. 

We note that an extrapolation technique similar to SSE has been used in Shell Model 
Monte Carlo calculations for over two decades 23|, [2^ . In that case, the extrapolation is 


performed by decomposing the Hamiltonian into “good sign” and “bad sign” parts, Hq and 
respectively. The calculations are then performed by multiplying the coefficients of 
by a parameter g and extrapolating from < 0 , where the simulations are free from sign 
oscillations, to the physical point g = 1- For SSE, the analysis in terms of “good” and “bad” 
signs is not the entire story. Most of the interactions can be divided into two groups which 
are “sign free” by themselves, such that a large portion of the sign oscillations is due to 
interference between the different underlying symmetries of the two groups of interactions. 
Since this effect is quadratic in the interfering interaction coefficients, the growth of the sign 
problem is more gradual. We therefore expect to be able to extrapolate from values not so 
far away from the physical point d^, = 1 . 

The LO auxiliary-held transfer matrix contains pion interactions with the matrix 

structure (TgTj acting on single nucleon states, and smeared contact interactions with matrix 
structures 1, itj and iugTj. Since 


0-2 G U:^[agT2\, U^[iag], ^^[icT^rg], 

0-2X3 e ^*[0-5X4], Wjixg], U^liagT^], 


(13) 

(14) 


we note that if we have an initial state with an even number of neutrons paired into spin- 
singlets and an even (but in general different) number of protons paired into spin-singlets. 
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then there will exist an antisymmetric representation for both (T 2 and (T 2 T 3 on the single 
nucleon states. The only interaction matrix structures in M^oaux included in both 
Eqs. (ITT]) and (1TT|) are and cr^Tg. Hence, the sign oscillations in det(M') will be produced 
by and (XgTg, and the interference between the two sets of interactions in Eqs. (ITT]) 
and (1141) . For initial states where the neutrons and protons cannot be paired into spin- 
singlets, there will be additional sign oscillations due to these unpaired nucleons. However, 
the number of such unpaired nucleons can often be kept to a minimum. 

III. RESULTS FOR CARBON-12 

We shall hrst discuss our results for the nucleus, as this provides us with a convenient 
test case for the SSE method. We extend the PMC calculation to larger values of the 
Euclidean projection time than would otherwise be possible, and verify how well these 
new results agree with earlier calculations. Thus, we shall work at hnite Euclidean time, 
extrapolate such data to —?• 1, after which the extrapolation ^ 00 is performed. For 
the case of the ®He nucleus, we shall consider the opposite order of limits and discuss in 
which cases either method is preferable. 

Our strategy for data analysis is to perform a global £t to all PMC data obtained for 
different values of dj^ and the coupling constant C 4 of the SU(4) symmetric Hamiltonian. 
The ansatz is 

n 

X(d,„ n) = .Y„ + .¥„"”<'■>(! - d,) + ^ A'f(15) 

i=i 

where the parameters Xq and are determined by a weighted least-squares fit. The 

superscript “SU(4)” implies that such parameters depend on C 4 . On the other hand, Xq 
represents the extrapolated value of the observable at = 1 , is therefore by dehnition 
independent of C^, and provides an important constraint for the least-squares fit. The de¬ 
termination of Xq for different observables provides the physical information of the analysis, 
and is obtained by a simultaneous fit to multiple instances of the (unphysical) coupling 
constant C 4 . In general, the results for d^j 7 ^ 1 depend on the choice of C^. 

The “order parameter” n in Eq. (IT5l) is adjusted to produce a per degree of freedom 
~ 1. For most observables, we find that n = 1 is sufficient, with some notable exceptions 
that we shall describe later. For n > 4, the per degree of freedom quickly saturates 
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below 1 and our fits become over-constrained. When different values of n produce similar 
per degree of freedom, we choose the lowest order for which we perform the extrapolation 
1. Our htting function ([13]) is suggested by the expectation that the dependence of 
observables on should, according to perturbation theory, be linear around d^ ~ 0 and 
df^ ~ 1. The linearity around d^ = 1 is incorporated by the sine function in Eq. fllSp . 

In Figs. [T]|U we show PMC data and extrapolations in df^ for the leading order (LO) con¬ 
tribution, the next-to-leading order two-nucleon force (NLO), the electromagnetic and strong 
isospin breaking (EMIB) and the next-to-next-to-leading order three-nucleon force (3NF), 





dh dh 


FIG. 1. PMC data at = 12.5 for with the upper figures showing the LO energy and 
the lower figures the shift due to two-nucleon forces at NLO. The left column shows the entire 
range in d^, while the right column shows a close-up around the physical point d^ = 1. The 
red lines are a simultaneous fit to all the PMC data using Eq. m, and each line corresponds 
to the Hamiltonian ()lip with a different choice of the coupling constant C^. The result of the 
extrapolation d^ —?• 1 is given by the red points in the right column of plots, along with the 
statistical uncertainty of the extrapolation. The error bars of the individual data points at d/^ < 1 
represent the Monte Carlo uncertainties. The cyan error bands correspond to the 67% confidence 
levels of the extrapolations. 
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FIG. 2. PMC data at = 12.5 for with the upper figures showing the total contribution 
from electromagnetic and strong isospin breaking (EMIB) operators, and the lower figures showing 
that of the three-nucleon force (3NF) operators at NNLO. Notation and conventions are identical 
to Fig. [H 

for the ground state energy at = 12.5 and = 14.5. The PMC data points in these 
hgures are comprised of runs with different choices of in the underlying SU(4) symmetric 
Hamiltonian. These are = —3.2 x 10“^, —3.4 x 10“®, —3.8 x 10“®, —4.2 x 10“^ and 
—4.8 X 10“^ (in units of MeV“^), which can be distinguished as separate bands in the hg¬ 
ures. The uppermost band corresponds to = —3.2 x 10“® MeV“^, and the lowest band 
to (^4 = —4.8 X 10“^ MeV“^. With the exception of the NLO contribution, all hts were 
performed with n = 1. For the NLO contribution, a higher order £t of n = 3 was used, in 
order to avoid introducing a bias around ~ 1 due to the increasingly accurate data at 
small values of df^. The higher order £t function has sufficient freedom to fully account for 
the data at small df^. 

From our PMC results, we hnd that even for large Euclidean projection times N^, the 
uncertainties of our calculations are well under control for < 0.8, which shows that the sign 
problem is under control for such values of d^. However, as d^^ —)■ 1 the uncertainties clearly 
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FIG. 3. PMC data (LO and NLO contributions) at = 14.5 for For a full description of the 
data and analysis, see the caption of Fig. [TJ 


grow aggressively, such that the sign problem reaches full strength at = 1 . As previous 
work has demonstrated, PMC calculations for at = 1 become very impractical for 
> 10. How far in the SSE analysis can be carried out depends on how robust the 
extrapolation in df^ can be made. We shall hnd that our extrapolated results remain robust 
as long as PMC calculations can be carried out for df^ > 0.75, and that this range could 
possible be extended to smaller df^ by choosing C 4 such that the extent of the linear region 
around = 1 is maximized. 


For ^^C, we have performed PMC calculations using SSE for 
between = 9.0 and 14.5 for one of the trial states used in Ref. 


iluchdean projection times 
j 8 | . It should be noted that 

the results of Ref. js] correspond to d;^ = 1, and could thus only be extended to ~ 10 
before the sign problem became prohibitive. In Fig. [5], we extend the dataset of Ref. 8 | 
with the new SSE data, and in Fig. [ 6 ] we combine all of our new results for ^^C with those 
originally shown in Ref. j^. We have also repeated the extrapolation in Euclidean time 
with this updated data, with a comparison and consistency check given in Table [B We hnd 
that our updated extrapolation is in agreement with the previous results of Ref. [81], and in 
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FIG. 4. PMC data (EMIB and 3NF contributions) at = 14.5 for For a full description of 
the data and analysis, see the caption of Fig. [2j 


most cases with a slightly reduced uncertainty. While the reduction in uncertainty is small, 
this is due to the fact that our new dataset is quite limited in comparison to that used in 
Ref. [St]. Moreover, the main objective of our study of is to establish that the SSE data 
are consistent with those obtained at = 1, when such calculations are not prohibited by 
the sign problem. 

It is worthwhile to investigate the stability of our extrapolated results as successive data 
points near = 1 are removed. In Fig. [71 we show the results of such a study for the 
ground state energy of at LO. Each data point in Fig. [7| corresponds to an extrapolation 

^ 1, such that the placement of the data point on the horizontal axis indicates the “break 
point” at which data is excluded from the analysis. For example, for an extrapolated value 
placed on the horizontal axis at d^ = 0.8, all PMC data with d^ > 0.8 have been excluded. 
For guidance, we also show a constant £t of the extrapolated points with 90% conhdence 
bands. The left panel of Fig. [7] shows data at = 12.5, while the right panel corresponds 
to = 14.5. In both cases, the extrapolated points located on the horizontal axis at d^^ = 1 
are the ones shown in FigjOl 
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Nt N, Nt N, 


FIG. 5. Comparison of the new PMC data for ^^C from the SSE analysis (red filled squares) and 
previous calculations ^ for df^ = 1 (blue open squares). The notation for the various contributions 
to the ground state energy E ^2 coincides with that of Table [H The results correspond to a trial 
state with an SU(4) coupling of —7.0 x 10“® MeV“^, not to be confused with the SU(4) coupling 
for the SSE analysis. It should be noted that the exponential deterioration of the Monte Carlo 
error has been circumvented. Also, these data should not be interpreted in terms of a “plateau” as 
a function of N^. An analysis of the dependence on is given in Fig . O and a concise description 
of the Euclidean time extrapolation method can be found in Ref. [25]. 


TABLE I. Contributions to the ground state energy of ^^C after extrapolation to inhnite Eu¬ 
clidean projection time. The contributions from the improved leading order amplitude (LO), the 
two-nucleon force at next-to-leading order (NLO), the electromagnetic and strong isospin break¬ 
ing (EMIB) and the three-nucleon force at next-to-next-to-leading order (3NE) are shown sepa¬ 
rately. The left column shows the results using the PMC data for = 1 from Ref. Qj, while the 
right column shows the results when the SSE data from this work are included. 



Ref. [8] 

+ SSE 

LO 

-96.92(16) 

-96.85(14) 

NLO 

10.48(3) 

10.47(3) 

EMIB 

7.76(1) 

7.76(1) 

3NE 

-14.80(6) 

-14.56(4) 
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FIG. 6 . Updated extrapolation in Euclidean projection time for the ground state. The new 
results from this work (red filled squares) have been combined with data compiled at = 1 from 


Ref. 


. The different data sets correspond to trial states with different SU(4) couplings (in units 


of MeV not to be confused with the SU(4) coupling C 4 for the SSE analysis. The results of the 


old and new ana 
detail in Ref. 251 


yses are given in Table HI and the extrapolation in Euclidean time is discussed in 


A determination of the smallest value of df^ from which a reliable extrapolation —)■ 1 
can be expected is highly significant, as systems with more nucleons and unequal numbers 
of protons and neutrons will suffer from an increasingly severe sign problem. In order to 
compensate, this forces PMC calculations to be performed at successively smaller (we note 
that for the sign problem is already apparent in the larger uncertainties at df^ = 0.95). 
As more complex nuclear systems are studied, the extent in which the extrapolated results 
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FIG. 7. Stability of the ground state energy at LO at = 12.5 (left panel) and = 14.5 
(right panel) under extrapolation 1. A data point placed at a given value of denotes that 

only data with equal and smaller values of are included in the extrapolation. Thus, a data point 
at d^ = 0.95 only includes data with d^ < 0.95. The horizontal bands with 90% confidence levels 
and are provided for visual guidance only. 


can be ascertained from lower valnes of dy^ will become more of an issne. Stndying this 
behavior with onr data provides an initial rongh idea on the robnstness of the SSE 
method. As can be seen from Fig. [71 onr extrapolated resnlt is in good agreement for valnes 
of df^ as low as ~ 0.8. Extrapolations nsing data below this valne only become increasingly 
nnreliable. This snggests the range in applicability of the dial parameter is limited. We hnd 
a similar conclnsion within the PMC calcnlations of ®He, which we shall tnrn to next. 


IV. RESULTS FOR HELIUM-6 AND BERYLLIUM-6 


The sign problem in the A = 6 system with 2 protons and 4 nentrons (or vice versa) 
is somewhat more severe than for Hence, if calcnlations are performed entirely at 
dj^ = 1, the extrapolation to inhnite Enclidean time (while still feasible) has to be performed 
nsing data with a rather limited range in Ny. However, for < 1 this sitnation improves 
rapidly. For ®He and ®Be, we shall therefore approach the problem differently than for 
We perform the extrapolation in Enclidean time for each pair of and d^, with the 
extrapolation d^ —)■ 1 as the hnal step of the analysis, in terms of Eq. flT5]) . For each valne 
of U 4 and d^j, a constrained Enclidean time extrapolation is performed nsing a minimnm of 
three trial states. This “triangnlation” strategy is described in detail in Ref. 2^. Once the 
data have been extrapolated to inhnite Enclidean time, the resnlting data are then snbjected 
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to the constrained global fitting procedure, as described for In Figs. | 8 ] and [H we show the 
£t results for the LO, NLO, EMIB, and 3NF contributions to the ®He energy. The explicit 
results are summarized in Table mi We hnd that the PMC results and the extrapolations 
for y4 = 6 are quantitatively similar to those of the system. The extrapolations of our 
results to = 1 are in good agreement with direct calculations at = 1 (where the sign 
problem is maximal) and we also hnd no indication of a breakdown or inconsistencies in the 
Euclidean time extrapolations as d^j —)■ 1. 

TABLE II. Contributions to the ground state energy of the A = 6 system after extrapolation 
d^ —>• 1. The contributions from the improved leading order amplitude (LO), the two-nucleon 
force at next-to-leading order (NLO), the electromagnetic and strong isospin breaking (EMIB) and 
the three-nucleon force at next-to-next-to-leading order (3NF) are shown separately. The leftmost 
column shows the result of a direct calculation at d^ = 1 without extrapolation in d^, and the other 
columns give the extrapolated results when progressively more data is excluded in the vicinity of 
dh = 1 - 



dh = 1 

4 < 0.95 

< 0.90 

< 0.85 

d,^ < 0.75 

LO 


-27.49(7) 

-27.54(4) 

-27.56(6) 

-27.56(10) 

-27.34(46) 

NLO 


2.61(4) 

2.58(2) 

2.61(3) 

2.66(4) 

3.05(11) 

EMIB 

(^He) 

1 .021(6) 

1.014(3) 

1.014(5) 

1.012(9) 

1.04(4) 

EMIB 

(®Be) 

2.65(1) 

2 .66(1) 

2.67(2) 

2.68(3) 

2.68(14) 

3NF 


-3.77(3) 

-3.76(1) 

-3.77(2) 

-3.73(2) 

-3.69(10) 


For the A = 6 system, we have performed PMC calculations for hve values of C^, namely 
-3.6 X 10-^ -4.2 X 10-^ -4.4 x 10-^-4.8 x 10-^ and -5.4 x lO"® (in units of MeV-^). 
We have employed a much larger range in d^j, and moreover each PMC data point has an 
order of magnitude better statistics compared with the data for ^^C. In all hgures related 
to A = 6 , the uppermost band of data corresponds to C 4 = —3.6 x 10“® MeV“^ and the 
lowest one to = —5.4 x 10“® MeV“^. As the data at low d^ exhibit increased curvature, 
the order n of our fit function has been taken to be n = 3 in contrast to the ^^C data, 
for which n = 1 was sufhcient in most cases. For ®Be, the LO, NLO, and 3NF results are 
identical to those shown in Figs. [8]and[9l However, the electromagnetic part of the EMIB 
contribution changes due to the different numbers of protons and neutrons. In Fig. [10], we 
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FIG. 8. PMC results for the ground state energy of ®He at LO (upper panels) and the contribution 
from two-nucleon interactions at NLO (lower panels). The left panels show the full range between 
= 0 and = 1, while the right panels show a close-up near d^ = 1. Each data point has 
been individually extrapolated to infinite Euclidean projection time before the SSE extrapolation 
d^ —)• 1. The results of the SSE extrapolation are shown by red data points in the right panels. 

show our results for the EMIB contribution to the energy of the ®Be ground state. 

As for we have also studied the stability of the extrapolated results for A = 6 as 
successive data points in the vicinity of = 1 are omitted. These results are summarized 
in Table HIl and the behavior of the extrapolated LO energy is illustrated in Fig. [TT] Our 
hndings suggest that a satisfactory extrapolation —)■ 1 can be obtained as long as PMC 
calculations can be performed for values no smaller than d^ ~ 0.80, at least for the present 
range of coupling constants employed in the analysis. We note that a larger range in 
may allow smaller values of d^ to serve as a useful starting point for the extrapolation, 
especially if the linear region around d;^ = 1 is thereby expanded. These hndings appear 
consistent with our conclusions for 

In contrast to performing extrapolations from SSE data below a given value of d/^, we can 
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FIG. 9. PMC results for the electromagnetic and strong isospin-breaking (EMIB, upper panels) 
and NNLO three-nucleon force (3NF, lower panels) contributions to the ground state energy of 
®He. The left panels show the full range between = 0 and = 1, while the right panels show 
a close-up near d^ = 1. Each data point has been individually extrapolated to infinite Euclidean 
projection time. Notation and conventions are as for Eig. [8l 



EIG. 10. PMC results for the electromagnetic and strong isospin-breaking (EMIB) component of 
®Be. The left panel shows the full range between d^ = 0 and d^ = 1, while the right panel shows 
a close-up near d^ = 1. Each data point has been individually extrapolated to infinite Euclidean 
projection time. Notation and conventions are as for Eig. [8j 
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FIG. 11. Upper panel: Stability of the ®He ground state energy at LO under extrapolation 1. 

The data point at = 1 represents a combined analysis including all the SSE data for d^ < 1 as 
well as the direct calculation at d^ = 1 without SSE. The data point at d^ = 0.95 only includes 
SSE data with d^ < 0.95 etc. Lower panels: Ground state energy of ®He at LO for d^ —>■ 1, using 
SSE data with d^ > 0.8 (left panel) and d^ > 0.7 (right panel). These data are the result of a 
linear extrapolation —>■ 1, such that each entry corresponds to one of ten possible datasets, as 
discussed in the main text. The horizontal band with 90% confidence level is provided for visual 
guidance only. 


instead perform linear fits to SSE data above a given df^. Snch an analysis can demonstrate 
in what range of a linear description remains valid in the vicinity of = 1. For this 
pnrpose, we perform linear fits for all possible combinations of three G 4 snbsets (of the five 
total). This provides ten different possible datasets. We then choose a valne of d^ below 
which all data is discarded, and perform a linear extrapolation to the remaining data points. 
In Eg. [TT], we show the extrapolated LO resnlts for > 0.7 and d^ > 0.8. We also show the 
average of the hts and corresponding 90% conhdence bands. For the analysis with dy^ > 0.8, 
onr average is consistent with the complete LO extrapolation shown in Table [TTl while for 
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FIG. 12. The mean value of the exponential of the complex phase (e*®) of det(iVT) as a function of 
Euclidean time step and SSE parameter for = —4.8 x 10“® MeV“^ in the A = 6 system. 
The interpolating lines are intended as a guide to the eye. 

the analysis with > 0.7, a clear systematical error is found. This shows (at least for ®He) 
that a linear description is valid for dj^ > 0.8. 

As an explicit demonstration of the amelioration of the sign problem using the SSE 
method, we show in Fig. [12] for A = 6 the dependence of the mean value of the exponential 
of the complex phase (e*®) of the PMC calculation on the number of Euclidean time steps 
and the SSE parameter df^ for = —4.8 x 10“^ MeV“^. For the case of df^ = 1, the 
mean value quickly approaches zero as is increased, indicating that the sign problem is 
becoming severe at rather modest Euclidean projection times. For decreasing values of d^, 
the effect of the sign problem is successively diminished, allowing the PMC method to be 
extended to signihcantly larger values of N^. 


V. DISCUSSION 

The SSE method introduced here is inspired by the existence of an SU(4) symmetric 
Hamiltonian which provides a reasonably accurate description of the physics of the full Chiral 
EFT Hamiltonian. This has already proven useful in earlier work, as it greatly facilitates 
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finding an accurate initial wave function which minimizes the extent of Euclidean time 
projection necessary with the full Hamiltonian. We have here illustrated how this concept 
can be taken one step further, by studying a weighted sum of physical and SU(4) symmetric 
Hamiltonians. In this way, the sign problem could be arbitrarily ameliorated, at the price 
of introducing an extrapolation in a control parameter —)■ 1. In practice, this means that 
the SSE method is only useful as long as the extrapolation errors can be kept under control. 
Naturally, performing simulations at a range of values of has the potential to multiply the 
required CPU time by a large factor. However, we have found that we are able to avoid an 
exponential increase in computation time as a function of Euclidean projection time, as long 
as we are able to perform simulations for > 0.75, as the accuracy of extrapolation then 
remains comparable with the statistical errors of typical simulations at = 1. We have also 
explored the freedom in the choice of the SU(4) symmetric Hamiltonian, which clearly plays 
no role at d^^ = 1, but which in general gives different results for df^ ^ 1. We have therefore 
made use of a “triangulation” method to improve the accuracy of the extrapolation d^j —)■ 1. 

An important consideration is whether a continuous shift from an SU(4) symmetric 
Hamiltonian to the full Chiral EFT Hamiltonian can be effected without inducing non¬ 
trivial changes in the spectrum. For instance, the appearance of a level crossing at a critical 
value of df^ would clearly limit the applicability of SSE. We note that such level crossings 
as a function of df^ would be quite rare for low-lying nuclear bound states, and furthermore 
we have the freedom to choose (7*4 to avoid such level crossings. But if a level crossing were 
to occur, there would be some subtleties in obtaining accurate and converged results. For 
such cases, it would be preferable to take the d;^ —)■ 1 limit first, followed by extrapolation 
in Euclidean time. An even better solution would entail solving a coupled-channel problem 
using multiple initial states. This makes it possible to disentangle one or more nearly de¬ 
generate states with the same quantum numbers. So far, we have concentrated on systems 
where such degeneracies are not expected, and where PMC is still possible (though difficult) 
without the SSE method. We have also not yet explored different ’’extrapolation Hamilto¬ 
nians” , where the sign oscillations are minimized while retaining as much as possible of the 
full Chiral EFT structure. 

In this first report on the SSE method, we have presented extensive results for ®He and 
^^C, for a wide range of df^ and SU(4) symmetric Hamiltonians. For the case of ®He, we have 
first extrapolated all results to infinite Euclidean time, followed by an extrapolation d^^ —)■ 1. 
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For the case of we have performed the extrapolation —)• 1 directly for data at hnite 
Euclidean projection time. In both cases, we hnd similar behavior as a function of df^ and 
excellent linearity in the vicinity of = 1. Which ordering of the limits is preferable depends 
on how severe the sign problem is for a given system at = 1. In cases where the sign 
problem is severe, the method of hrst taking the limit d^ —)■ 1 followed by an extrapolation 
in Euclidean time is preferable. This is because the Euclidean time extrapolation may be 
ambiguous if only very short Euclidean times are accessible with PMC. Our results provide 
the hrst determination of the ground state energy of ®He using lattice Chiral EFT, and for 
^^C we hnd that our new results at larger Euclidean projection time agree very well with 
previous inhnite-time extrapolations. 


We have presented here the hrst studies of the binding energy of ®He within lattice Chiral 
EFT. While the results are very encouraging for the feasibility of calculations for neutron- 
rich systems, we also hnd that ®He appears underbound by ~ 1 MeV at NNLO. Nevertheless, 
since the present calculations are performed in an L = 6 box, the possibility remains that 


hnite volume ehects could signihcantly improve on t 


F-wave states are underbound in a hnite volume 


le current situation. In particular, since 


26 


27| while F-wave states are overbound. 


calculations in larger boxes (or extrapolations to L = cxd) would shed more light on this 
situation. 
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